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HIGH-LEVEL  SPACECRAFT  CHARGING  AT  GEOSYNCHRONOUS  ALTITUDES: 

A  STATISTICAL  STUDY 
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Maurice  Tautz 
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Abstract 


We  present  the  results  of  a  statistical  study  on  high-level  spacecraft  charging  at 
geosynchronous  altitudes.  Below  the  critical  temperature  T*  for  a  surface  material,  no  spacecraft 
charging  occurs.  The  spacecraft  charging  potential  data  are  obtained  from  the  ion  line  of  the  ion 
energy  spectrum.  If  the  ion  line  can  not  be  clearly  identified,  the  data  point  is  flagged.  We  do  not 
use  flagged  data  in  our  analysis.  Since  T*  depends  on  the  surface  material  and  since  each 
satellite  has  its  own  surface  material  or  materials,  each  satellite  is  expected  to  have  its  own 
critical  temperature.  The  coordinated  space  environmental  parameter  data  of  the  Los  Alamos 
National  Laboratory  (LANL)  geosynchronous  satellites  include  spacecraft  charging  data 
measured  on  several  geosynchronous  satellites  in  eclipses  and  in  sunlight  over  years.  We  have 
obtained  statistical  results  of  T*  for  each  satellite  studied  and  found  that,  beyond  T*,  the  high- 
level  spacecraft  potential  increases  almost  linearly  with  the  ambient  electron  temperature. 
Amazingly,  the  critical  temperature  in  sunlight  remains  the  same  as  in  eclipse,  agreeing  with  the 
monopole-dipole  differential  charging  model.  This  work  offers  a  useful  method  not  only  for 
predicting  the  onset  of  spacecraft  charging  in  eclipse  and  in  sunlight  but  also  for  predicting  high- 
level  spacecraft  charging  potential  with  reasonable  accuracy  at  any  given  ambient  electron 
temperature  in  the  geosynchronous  environment. 

Introduction 


Spacecraft  charging  can  detrimentally  affect  electrical  operations  on  space  systems.  Most 
communication  and  surveillance  spacecraft  are  at  geosynchronous  altitudes  and  many  more  such 
spacecraft  will  be  deployed  in  the  new  millennium.  The  plasma  density  in  the  geosynchronous 
environment  varies  from  over  100  cm'3  to  0.1  cm'3  and  the  energy  varies  from  a  few  eV  to  tens  of 
keV  depending  on  local  time  and  geomagnetic  conditions.  Spacecraft  surface  charging  occurs  at 
high  plasma  energies.  While  surface  material  properties  and  spacecraft  geometry  are  defined  by 
spacecraft  design,  spacecraft  charging  is  controlled  by  the  dynamic  plasma  condition  which 
varies  in  time.  It  is  important  to  identify  the  most  reliable  space  environment  parameters  for 
predicting  spacecraft  charging. 

We  have  identified  the  most  reliable  space  environment  parameter  for  predicting  spacecraft 
charging  is  the  plasma  electron  temperature.  Historically,  Rubin  et  al.  [1980],  first  reported, 
though  with  a  few  data  points  only,  the  linear  dependence  of  spacecraft  potential  as  a  function  of 
the  ambient  electron  temperature.  Their  graph  of  spacecraft  potential  versus  electron 
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temperature  showed  an  intercept  at  a  finite  temperature.  Laframboise  et  al.  [1982],  Lai  et  al. 
[1982,  1983],  and  Prokopenko  and  Laframboise  [1983]  put  forth  the  theory  of  the  critical 
temperature.  Below  the  critical  temperature,  spacecraft  charging  does  not  occur;  above  it, 
spacecraft  charging  occurs.  The  theory  is  now  included  in  standard  textbooks  on  spacecraft- 
plasma  interactions  [Hastings  and  Garrett,  1996]. 

For  many  years,  there  was  no  systematic  observation  of  spacecraft  charging  together  with  the 
coordinated  space  environment  parameters  and  therefore  no  way  to  validate  the  theory. 

Recently,  the  spacecraft  charging  data  obtained  on  the  Los  Alamos  National  Laboratory  (LANL) 
geosynchronous  satellites  LANL-89-046,  LANL-90-95,  LANL-91-80,  LANL-94-084  and 
LANL-97A  have  become  available.  A  recent  study  [Lai  and  Della-Rose,  2001]  using  four  weeks 
of  spacecraft  charging  data  obtained  on  LANL-94-084  has  uncovered  new  evidence  for  the 
existence  of  critical  temperature  of  the  space  plasma  electrons  for  a  given  spacecraft  surface 
material.  Below  the  critical  temperature,  little  spacecraft  charging  occurs,  while  above  the 
critical  temperature,  the  spacecraft  potential  increases  almost  linearly.  The  evidence  that  a 
critical  temperature  exists  is  abundant  and  was  seen  in  every  charging  event  on  the  LANL 
satellite. 

In  order  to  establish  the  critical  temperature  theory  beyond  a  doubt,  it  is  necessary  to  study 
more  data  for  validation.  We  report  on  our  results  in  this  paper. 

Theory 

In  the  geosynchronous  environment,  the  ambient  electron  flux  exceeds  that  of  the  ambient 
ions  by  two  orders  of  magnitude,  because  the  electrons  are  lighter  and  faster.  Measurements  on 
the  LANL-94-084  satellite  confirmed  the  flux  difference  [Lai  and  Della-Rose,  2001]. 

Suppose  a  satellite  is  initially  uncharged  and  the  ambient  electron  temperature  is  increasing. 
Eventually,  the  electron  flux  increases  to  a  level  for  charging  to  occur.  At  the  threshold,  the 
current  balance  is  between  the  incoming  ambient  electrons  and  the  outgoing  secondary  and 
backscattered  electrons. 

f  dEEf(E)  =  f  dEEf(E)[S(E)  +  n(E)]  (1) 

*(E)  and  0(E)  are  the  coefficients  of  secondary  electron  emission  and  backscattered  electron 
emission  respectively.  For  a  Maxwellian  space  plasma,  the  distribution  function /(E)  is  of  the 
form: 

f(E)  =  n(m/2ttkT)'nQxp(-E/kT)  (2) 

Substituting  eq(2)  into  eq(l),  one  readily  obtains  two  theorems,  (I)  and  (II):  - 

(I)  Since  the  density  n  is  multiplicative,  it  cancels  out  on  both  sides.  Therefore,  the 
threshold  condition  is  independent  of  the  plasma  density  n. 

(II)  The  solution  T*  to  eq(l)  is  the  critical  temperature  for  the  onset  of  spacecraft 
charging. 

To  solve  eq(l),  one  needs  to  input  the  functions  *(E)  and  0(E).  If  the  ambient  electrons  are 
coming  in  at  various  angles  2,  one  needs  to  use  angle  dependent  functions  of  *(E,2)  and  0(E,2). 


Including  the  angles,  the  algebra  becomes  more  complicated.  The  results  are  given  in  Table  1 . 


Table  1.  Critical  Temperature  (keV) 


MATERIAL 

ISOTROPIC 

NORMAL 

Mg 

0.4 

— 

A1 

0.6 

— 

Kapton 

0.8 

0.5 

A1  Oxide 

2.0 

1.2 

Teflon  © 

2.1 

1.4 

Cu-Be 

2.1 

1.4 

Glass 

2.2 

1.4 

Si02 

2.6 

1.7 

Silver 

2.7 

1.2 

Mg  Oxide 

3.6 

2.5 

Indium  Oxide 

3.6 

2.0 

Gold 

4.9 

2.9 

Cu-Be  (Activated) 

5.3 

3.7 

MgF2 

10.9 

7.8 

Spacecraft  Potential 


For  the  onset  of  spacecraft  charging,  we  have  neglected  the  ions  because  the  ion  current  is 
about  two  orders  of  magnitude  smaller  than  that  of  the  electrons.  However,  to  determine  the 
resulting  spacecraft  equilibrium  (negative)  potential,  one  needs  to  include  the  ions  that  are 
attracted.  It  is  often  a  good  approximation  to  describe  the  electron  repulsion  by  means  of  the 
Boltzmann  factor  (the  exponential  term)  and  the  ion  attraction  by  means  of  the  factor  of  Mott- 
Smith  and  Langmuir  (the  term  in  parenthesis  following  the  ion  current  7,  (0))  in  the  following 
current-balance  equation: 


ZW)  =  /, (*)-/.(#)  =  /,(' 0) 


1-- 


kT 


-4(0)  exp 


i  J 


kT 


=  0 


(3) 


where  Ie  (A)  and  7j  (A)  are  the  currents  of  the  repelled  and  attracted  species,  respectively, 
collected  by  the  spacecraft  at  potential  N.  Including  secondary  and  backscattered  electrons,  the 
electron  current  7e(N)  in  eq(3)  is  replaced  by  the  net  electron  current ,  which  is  of  the  form: 


lM,Ts)  =  W)[l-<d+t}>]  (4) 

where 

f dEEf(E,Te)[S(E)  +  rtE )] 

<S  +  tj>=— - - -  (5) 

£ dEEf(E,Te ) 

Using  eq(5),  the  current-balance  equation,  eq(3),  becomes 


f  ej) 

f 

4(0) 

1~  ~  4  (0)  [1-  <  £ + 77  >]  exp 
l  kTi) 

r«rj 

where  the  subscripts  e  and  i  label  electrons  and  ions,  ee  =  -e  and  ei  =  e,  and  e  is  the  elementary 
charge.  For  spacecraft  potential  lower  than  the  ambient  ion  temperature,  that  is  qN«  kT  j, 


which  is  usually  valid  initially,  one 


e<j>  = 


TIT 


■kTe  log 


/,(0) 


+  ... 


T,/Te+\  e  ~Ie(0) 

which  gives  a  straight  line  of  eNas 
needs  to  be  included. 


expands  eq(3.6)  in  a  Taylor  series  of  qNI  kT ,  and  obtains 

(7) 

a  function  of  kT.  For  higher  values  of  eN,  the  quadratic  term 


Observations  of  Spacecraft  Charging  in  Eclipse 

The  ion  to  electron  temperature  ratio  Tj/Te  obtained  on  LANL1994-084  is  mostly  between  1 
and  3  approximately.  If  the  ratio  is  almost  a  constant  in  eq(7),  the  potential  //plotted  as  a 
function  of  kT e  would  be  almost  a  straight  line.  It  is  not  really  a  straight  line,  because  the  log 
term  is  also  a  function  of  T  e. 


The  LANL  spacecraft  potential  data  are  deduced  from  the  shift  of  the  ion  distribution /i(E). 
Occasionally,  the  ion  shift  (ion  line)  is  not  clear.  Instead  of  discarding  those  points  which  have 
unclear  ion  shifts,  the  LANL  team  deduced  the  spacecraft  potential  for  those  points  in  the 
following  manner  [M.F.  Thompsen,  personal  communication.,  2002].  They  have  observed  an 
empirical  curve  of  spacecraft  potential  plotted  against  electron  temperature,  although  they  have 
made  no  mention  of  critical  temperature.  They  put  the  ‘unclear’  points  on  this  empirical  curve 
according  to  their  electron  temperature.  Such  data  are  flagged.  We  have  examined  the  flagged 
data  points  in  three  different  months  of  sample  data  [courtesy,  M.F.  Thompsen,  2002]  and  found 
that  the  flagged  data  are  all  on  the  sunlight  charging  curve  and  none  on  the  eclipse  charging 
curve.  Therefore  we  conclude  that  the  flagged  data  probably  do  not  affect  our  eclipse  charging 
study  at  all. 


Figure  1  Spacecraft  potential  and  electron 
temperature  measured  on  L AN L- 1 989-046 
during  the  eclipse  periods  of  Mar  13-28, 1993- 
2001. 


4 

ELECTRON  TEMPERATURE  (keV) 


Figure  2  Spacecraft  potential  and  electron 
temperature  measured  on  LANL- 90-095 
during  the  eclipse  periods  of  Mar  14-29, 
1993-2001. 


Figure  1  shows  a  plot  of  LANL  1989-046  spacecraft  potential  vs.  electron  temperature  for 
the  eclipse  periods  of  Mar  13  to  28,  1993  to2001.  Figure  2  shows  a  similar  plot  of  LANL  1990- 
095  spacecraft  charging  and  ambient  electron  temperature  data  taken  in  the  eclipse  periods  of 
Mar  14  to  29, 1993  to  2001.  We  have  also  obtained  similar  results  by  using  eclipse  data  for  the 
LANL- 1994-084,  LANL  1991-80  and  LANL-1997A  satellites.  Each  of  these  results  shows  that 
(1)  the  trend  of  each  plot  is  almost  a  straight  line.  (2)  there  exists  an  intercept  T*  (its  value  being 
1.5  to  2  keV  approximately),  (3)  a  quadratic  function  would  fit  better  than  a  straight  line,  (4)  the 
values  of  T*  are  lower  with  quadratic  fits,  but  are  of  the  same  order  of  magnitude. 


Figure  3.  The  spacecraft  potential  mirrors  with  the  electron  temperature  -  not  with  the 
electron  density  or  kp. 

The  most  well  known  storm  in  the  recent  solar  max  was  the  Bastille  Day  storm  of 2000. 
Figure  3  shows  the  magnetic  index  kp  going  up  to  the  maximum  possible  level  (k=9)  during  the 
storm.  The  parallel  and  perpendicular  electron  temperatures  jumped  to  high  values  at  the  arrival 
of  the  solar  disturbance.  The  electron  density  rose  later  but  not  simultaneously  with  the  electron 
temperature.  This  offers  a  good  opportunity  to  observe  the  spacecraft  potential  in  response  to 
these  two  space  parameters,  namely,  electron  temperature  and  electron  density.  Indeed,  the  data 
measured  on  LANL  1994-084  shows  that  the  spacecraft  potential  mirrored  the  electron 
temperature  but  not  with  the  electron  density  at  all.  It  had  been  previously  assumed  that  electron 
flux  and  electron  density  were  the  principal  factors. 

Spacecraft  Charging  in  Sunlight 

Photoelectron  emission  from  spacecraft  surfaces  often  dominates  over  all  other  ambient 
currents.  Therefore,  it  seems  impossible  for  spacecraft  to  charge  to  negative  potentials  in 


sunlight.  Yet,  the  LANL  satellites  often  charge  to  hundreds  or  thousands  of  negative  volts  in 
sunlight.  We  are  going  to  present  that  the  LANL  satellite  charging  data  agree  well  with  the 
monopole-dipole  model. 

In  the  following,  we  will  present  an  important  finding:  “The  level  of  spacecraft  charging  in 
sunlight  is  usually  about  1/3  of  that  in  eclipse”.  As  a  collorary,  we  have  another  important 
finding:  “The  onset  of  spacecraft  charging  occurs  at  the  same  critical  temperature  no  matter  the 
satellite  is  in  clipse  or  in  sunlight.” 

Photoelectron  Current  in  Sunlight 

Laboratory  measurements  show  that  typical  spacecraft  surface  materials  generate  in  space 
conditions  photoemission  currents  of  the  order  of  Jph  =  2x1  O'9  A/cm2  [Stannard,  et  al.,  1978].  For 
comparison,  the  average  electron  current  density  measured  on  the  SCATHA  (Spacecraft 
Charging  at  High  Altitudes)  satellite  was  <J>  =  0.1 15xl0'9  A/cm2  [Purvis  and  Garrett,  1984]. 
Thus,  the  photoemission  current  density  Jph  exceeds  the  average  ambient  electron  current  J  by  a 
factor  F  =  20. 

Iph~20<J>  (8) 


If  a  spacecraft  is  charged  to  a  negative  potential  in  steady  state,  the  potential  must  satisfy  the 
current  balance  equation.  For  example,  current  balance  in  the  Mott-Smith  and  Langmuir  model 
is  of  the  form: 


1-- 


k T 


=  Jph  (9) 


V  nt-i  J 

where  the  notations  are  as  usual,  qe=  -e,  q\  =  e,  and  D  is  negative.  If  /Ph  >  /e(0),  there  is  no 
solution  for  the  above  equation.  How,  then,  can  spacecraft  charging  to  negative  potentials 
occur? 


Monopole-Dipole  Potential  in  Sunlight 


Consider  a  negatively  charged  satellite  with  dielectric  surfaces  in  sunlight.  The  shadowed 
side  charges  to  a  higher  (negative)  potential  than  the  sunlit  side.  Taking  the  first  two  terms  of  an 
infinite  series  in  Legendre  function  of  a  general  potential  distribution,  one  obtains  a  monopole- 
dipole  potential  distribution  [Besse  and  Rubin,  1980]: 

(l0) 


where  A  is  the  dipole  strength  which  is  less  than  unity.  The  angle  2  =  0°  is  the  normal  sunlight 
direction  (Figure  4).  We  have  conducted  an  explicit  validation,  for  the  first  time,  of  the 
monopole-dipole  distribution  in  sunlight  charging. 


For  unit  radii  satellite,  the  radial  distance  /?=1  at  the  satellite  surface,  and  K  =  N(90°,  1) 
equals  the  monopole  potential.  The  potential  barrier  is  located  at  Rs,  where  the  potential  is  the 
maximum  for  2  =  0°. 


where  the  photoelectron  distribution  is  Maxwellian  and  has  a  typical  temperature  T  of  1.2  eV 
[Whipple,  1982].  Since  the  photoelectron  temperature  is  low,  very  small  barrier  height  will 
block  most  of  the  photoelectrons. 


Figure  4.  Typical  potential  contours  in  a  monopole-dipole  potential  distribution.  The 
contours  wrap  from  the  dark  side  to  the  sunlit  side. 


Figure  5.  Fraction  of  photoelectron  current  escaping  through  the  barrier.  A  small  barrier 
height  is  sufficient  to  block  most  of  the  photoelectron  current. 

In  particular,  for  high  level  charging,  the  ratio  B/K  of  the  barrier  potential  to  the  monopole 
potential  is  usually  nearly  zero.  Substituting  B/K  ~  0  in  eq(12)  gives  A  ~  'A.  As  a  result,  the 
ratio  of  the  sunlit  surface  potential  to  that  of  the  shaded  surface  is  given  by  (Figure  6): 

*(0M)  1  -A  1-1/2  1 

^(180°,  1)  l  +  A  1  +  1/2  3  1  ' 


Observations  of  Spacecraft  Charging  in  Sunlight 


We  have  analyzed  several  years  of  the  Los  Alamos  National  Laboratory  (LANL)  satellite 
charging  data  and  obtained  very  useful  statistical  results  on  the  behavior  of  charging  in  sunlight. 
The  results  agree  surprisingly  well  with  the  monopole-dipole  model. 

Figure  7  presents  a  compilation  of  charging  data 
obtained  on  the  LANL  1991-80  Satellite  for  the  period 
March  13-28,  1994-2001.  The  March  period  is  chosen 
because  of  the  following  reason.  The  geosynchronous 
satellite,  being  in  the  equatorial  plane,  undergoes  eclipse 
around  midnight  in  March  and  September.  Therefore, 
both  eclipse  charging  and  sunlight  charging  data  are 
available  for  comparison  during  the  March  period.  A 
statistical  curve  fit  shows  clearly  that  there  are  two  trends 
in  Figure  7.  The  ratio  of  the  potentials  in  eclipse  and  in 
sunlight  is  approximately  1/3.  For  example,  when  the 
sunlight  charging  potential  is  -1  kV  at  electron 
temperature  of  about  4.4  kV,  the  eclipse  charging 
potential  at  the  same  temperature  is  about  -3  kV  .  This 
result  is  a  strong  evidence  supporting  the  validity  of  the 
monopole-dipole  model. 


Figure  6.  Ratio  of  potentials  of 
sunlit  and  shadowed 
surfaces  as  a  function  of 
barrier  potential. 


Since  1/3  of  a  finite  number  is  finite  while  1/3  of  zero  is  obviously  zero,  the  onset  of 
spacecraft  charging  occurs  at  the  same  critical  temperature  of  the  space  plasma  electrons.  Below 
the  critical  temperature,  there  is  no  charging.  Above  it,  the  charging  level  increases  with  the 
temperature.  This  new  and  important  finding  is  clearly  seen  in  Figure  7. 

To  add  more  credence  to  the  above  findings,  we  present  another  statistical  graph  (Figure  8) 
showing  the  charging  data  of  a  different  satellite,  viz.,  LANL-1997A,  during  the  periods 
September  14-28, 1997-2001.  Again,  the  period  September  14-28  is  chosen  because  both  eclipse 
charging  and  sunlight  charging  events  are  available  for  comparison.  In  Figure  8,  the  ratio  of 
potentials  in  sunlight  to  that  in  eclipse  is  again  approximately  1/3.  For  example,  the  sunlight 
charging  potential  is  about  -1  kV  at  the  temperature  of  about  4.4  kV;  that  of  the  elipse  charging 
at  the  same  temperature  is  about  -3  kV.  The  onsets  of  eclipse  charging  and  sunlight  charging 
occur  at  approximately  the  same  value  of  critical  temperature.  Similar  results  have  been 
obtained  for  the  other  LANL  satellites.  We  conclude  that  these  results  strongly  support  the 
validity  of  the  monopole-dipole  model. 


ELECTRON  TEMPERATURE  (kaV) 

Figure  8.  Charging  voltage  in  sunlight  is  1/3  that 
in  darkness. 


Summary  and  Discussion 


The  geosynchronous  altitudes  are  the  most  important  region  for  spacecraft  charging.  Most 
communication  satellite  are  there,  while  the  ambient  plasma  varies  in  temperature  and  density 
depending  on  the  local  time  and  space  weather.  A  most  important  problem  is  to  seek  the  most 
reliable  space  environment  indicator  of  the  onset  of  spacecraft  charging.  By  studying  the  LANL 
geosynchronous  satellite  data,  we  have  found  that  spacecraft  charging  depends  strongly  on  the 
ambient  plasma  electron  temperature,  and  not  as  much  on  the  other  space  environment 
parameters  such  as  electron  density,  ion  temperature,  ion  density,  and  even  high  kp.  The  exact 
level  of  charging  depends  on  the  spacecraft  geometry,  configuration  of  surfaces,  etc.  The 
average  potential  is,  however,  almost  a  straight  line  as  a  function  of  temperature.  The  line 
deviates  from  being  absolutely  straight,  because  there  is  a  quadratic  term  in  the  Taylor  series  of 
potential  as  a  function  of  temperature  and,  furthermore,  the  plasma  distribution  may  resemble  a 
kappa  distribution  at  high  temperatures. 

It  has  long  been  known  from  laboratory  measurements  that  the  photoemission  flux  from 
sunlit  surfaces  should  exceed  the  ambient  plasma  electron  flux  at  geosynchronous  altitudes. 

How  then  can  satellites  charge  to  negative  voltages?  To  solve  this  problem,  one  realizes  that 
there  is  always  a  shadowed  surface  when  a  satellite  is  in  sunlight.  High  negative  charging  can 
occur  on  the  shadowed  surface  when  the  plasma  electron  temperature  is  hot. 

A  bootstrap  mechanism  then  occurs  as  follows.  The  negative  potential  contours  extend  from 
the  high  voltage  (dark)  side  to  the  low  voltage  (sunlit)  side,  thereby  trapping  the  photoelectrons 
with  a  potential  barrier.  Since  photoelectrons  have  a  few  eV  only  in  energy,  a  small  potential 
barrier  is  sufficient  to  trap  most  of  the  photoelectrons.  With  photoemission  partly  suppressed, 
the  sunlit  side  can  charge  to  high  negative  potentials.  The  results  of  the  monopole-dipole  model 
predicts  that  the  ratio  between  the  eclipse  charging  potential  to  the  sunlit  charging  potential 
under  the  same  space  plasma  conditions  is  approximately  1/3. 

Amazingly,  we  have  found  abundant  evidence  that  this  ratio  is  valid  statistically  by  studying 
years  of  charging  data  obtained  on  LANL  satellites, .  Since  R  of  a  finite  number  is  finite  and  R 
of  a  zero  is  zero,  the  critical  temperature  for  the  onset  of  spacecraft  charging  is  therefore 
unchanged  -  it  is  the  same  whether  in  eclipse  or  in  sunlight.  We  believe  that  this  is  a  significant 
finding  for  developing  a  reliable  indicator  in  the  future  for  predicting  spacecraft  charging  in 
changing  space  weather. 
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